p 


AD- A  1 34  458 


UNCLASSIFIED 


THE  EVOLUTION  OF  SPONTANEOUS  AND  COHERENT  RADIATION  IN 
THE  FREE  ELECTRON  LASER  OSCILLATOR ( U )  NAVAL  RESEARCH 
LAB  WASHINGTON  DC  P  SPRANGLE  ET  AL-  16  AUG  83 
NRL -MR -5 1 10  '  F/G  20/5 


1/1 


NL 


m  FILE  COPY  /4/JV  ys'P 


r 


<3> 

NRL  Memorandum  Report  Slid 

The  Evolution  of  Spontaneous  and  Coherent 
Radiation  in  the  Free  Electron  Laser  Oscillator 

P.  Sprangle,  C.  M.  Tang  and  I.  Bernstein* 

Plasma  Theory  Branch 
Plasma  Physics  Division 

*Yale  University 
New  Haven,  CT 


August  16,  1983 


This  work  was  supported  by  DARPA  under  contract  No.  3817. 


NAVAL  RESEARCH  LABORATORY 
Waskinctoo,  D.C. 


DTIU 


^LECTE 

NOV  0  7 1983 


Approved  for  pubfic  rekese;  dkoibutioa  unlimited. 


88  11  04  050 


SECURITY  C  L  ASSl  F  »  C  A  T  i  ON  Of  7  *»  i  S  f*  A  \>  f  •>>«#»  /)«(•  (m«i«j 


REPORT  DOCUMENTATION  PAGE  j 

)  REPORT  NUMBER  2  GOVT  ACaSSlON  l»Oi  1  C  iP.  1  N  ▼  ’  S  C  A  T  0-  nuMBE  R 

NRL  Memorandum  Report  5110  W  S'p 

4  T|TuC  (and  Subtitle)  ^ 

THE  EVOLUTION  OF  SPONTANEOUS  AND 
COHERENT  RADIATION  IN  THE  FREE 

ELECTRON  LASER  OSCILLATOR 

1  TrPE  of  REPORT  »  Pf»,(JD  COVERED 

Interim  report  on  a  continuing 
NHL  problem. 

6  PEA*  O pm i  s  G  ORG  REPORT  NuMBER 

7  AUTHOR/ •> 

P.  Sprangle,  C.M.  Tang  and  1.  Bernstein* 

- 

i  contract  or  grant  number/*) 

9  PERFORMING  ORGANIZATION  NAME  AnO  ADDRESS 

Naval  Research  Laboratory 

Washington,  DC  20375 

io  program  element  project  task 

AREA  6  WORK  UNlT  NUMBERS 

62301 E;  47-0867-0-3 

n  CONTROLLING  Of  fCE  NAME  »NO  AOORESS 

Defense  Advanced  Research  Projects  Agency 

Arlington,  VA  22209 

12  REPORT  Date 

August  16, 1983 

1)  NUMBER  OF  PAGES 

39 

14  MONITORING  AGENCY  nAME  A  ADDRESS/!/  different  horn  Controlling  Ottlce) 

IS  SECURITY  CLASS  (of  thle  report) 

UNCLASSIFIED 

ISa.  OECLASSlFl  CATION  7  DOWN  GRADING 
SCnEOULE 

16  DISTRIBUTION  STATEMENT  (ol  thl  m  Report, 

Approved  for  public  release;  distribution  unlimited. 

17  DISTRIBUTION  STATEMENT  /o/  /A*  abstract  entered  /n  Block  50.  It  dltleranl  from  Report) 


18  supplementary  notes 

♦Present  address:  Yale  University,  New  Haven,  CT 

This  work  was  supported  by  DARPA  under  contract  No.  3817. 

19.  KEY  WOROS  (Conf/nu#  on  raveree  aid*  //  necaeemry  and  identity  by  block  number) 

Free  electron  laser  Incoherent  radiation 

Spontaneous  radiation  Oscillator  start-up 


20.  ABSTRACT  (Continue  on  ravaraa  a/da  It  nacaaaary  and  Identity  by  block  number) 

- 

An  analysis  of  the  free  electron  laser  (FEL)  oscillator  start-up  problem  in  the 
linear  regime  is  presented.  The  model  is  spatially  one  dimensional,  though  many 
important  three  dimensional  effects  are  included  heuristically.  The  electron  beam 
consists  of  pulses  of  arbitrary  shape  separated  by  approximately  twice  the  radiation 
transit  time.  The  small  gain /pass  approximation  is  employed  in  deriving  an  energy 

(Continues) 

DO  I  JAn"?!  1^73  COITION  OF  I  NOV  f  IS  OMOLETI! 


S/N  0  10  2*  0  14*  660  I  - _ _ _ ___ 

SECURITY  CLASSIFICATION  OF  THIS  PACE  fBfcan  Data  *n  fared. 


1 


JO  ABSTRACT  (Continued) 


^  rate  equation,  which  described  the  evolution  of  the  radiation  pulses  within  the 
resonator.  The  wiggler  field  is  assumed  to  occupy  a  portion  of  the  finite  Q 
resonator.  In  the  energy  rate  equation,  the  spontaneous  (incoherent)  radiation 
term  is  represented  by  a  source  matrix,  while  the  stimulated  (coherent)  radiation 
term  is  represented  by  a  gain  matrix.  The  effect  of  small  variations  in  the 
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suggests  possible  methods  which  could  substantially  shorten  the  start-up  times  in 
FEL  oscillators.  Finally,  our  results  are  compared  with  the  FEL  oscillator 
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THE  EVOLUTION  OF  SPONTANEOUS  AND  COHERENT 
RADIATION  IN  THE  FREE  ELECTRON  LASER  OSCILLATOR 

I .  Introduction 

A  number  of  successful  free  electron  laser  oscillator  experiments  have  been 

reported.^  ^  Simple  considerations  concerning  the  spontaneous  radiation  level 

a 

indicated  start-up  times  much  shorter  than  those  observed.  Since  a  number  of 

experiments  utilizing  shorter  electron  beam  macropulses  are  being  constructed  or 

planned,  thus,  there  is  concern  that  these  forthcoming  experiments  may  be  unable 

to  reach  saturation.  A  quantitative  understanding  of  the  growth  of  coherent 

stimulated  radiation  from  incoherent  spontaneous  emission  is  thus  highly 

desirable.  Published  papers  on  the  FF.L  oscillator  had  either  neglected  the 

s_20 

spontaneous  radiation  ,  or  had  treated  them  separately  from  the  stimulated 
1 3 

radiation.  Here  we  outline  a  classical  theory  of  the  spontaneous  start-up  of 
the  FEL  oscillator  in  the  cold,  small  signal  regime.  Our  model  is  spatially  one 
dimensional  and  therefore  lacks  many  important  features  such  as  transverse 
gradients  associated  with  the  radiation  and  electron  beam  and  diffraction 
effects.  In  a  one-dimensional  model  these  effects  can  only  be  incorporated  in  an 

approximate  way  by  means  of  filling  factors. 

5~*28 

Theories  of  the  free  electron  laser  (FEL)  have  proceeded  from  a  continuum 

description  of  the  electron  dynamics,  either  fluid  equations  or  the  Vlasov 
equation.  For  a  proper  description  of  the  start-up  of  an  FEL  oscillator  one  must 
take  into  account  the  fact  that  the  electrons  are  discrete  and  initially 
uncorrelated,  since  it  is  the  acceleration  radiation  of  individual  electrons  in 
the  wiggler  that  provides  the  initial  fields.  These  initial  fields  are  then 
amplified  by  the  collective  gain  mechanism  associated  with  the  continuum 
description.  This  Initial  radiation,  however,  is  effectively  incoherent  in  a 

device  in  which  the  electron  density  is  small  and  the  electrons  are  randomly 
Manuscript  approved  May  23,  1983. 
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distributed.  Thus  a  statistical  theory  is  required  which  is  couched  in  terms  of 
objects  bilinear  in  the  fluctuating  quantities  so  that  ensemble  averages  are  non¬ 
zero,  even  when  the  ensemble  average  fluctuating  current  density  is  zero. 

The  theory  described  in  this  work  is  one  dimensional  in  space  and  treats  the 
electrons  as  governed  by  the  relativistic  equations  of  motion,  and  the 
electromagnetic  fields  as  governed  by  Maxwell's  equations.  This  is  valid  whenever 
the  root  mean  square  fluctuation  6N  in  the  number  of  photons  in  the  resonator  is 
small  compared  with  che  mean  number  of  photons  N.  Certainly  this  is  not  true 
initially,  and  in  principle,  one  should  treat  the  problem  initially  by  quantum 

mechanics.  Failure  to  do  so  implies  an  uncertainty  in  the  initial  phases  of  the 

-  1/2 

start  up.  Since  one  expects  6N  ~  (N)  if  the  electrons  are  randomly 

distributed,  the  duration  of  the  quantum  regime  will  be  short  if  classical  theory 

predicts  for  times  short  compared  to  that  for  saturation  that  the  photon  density 
3  2  2  2 

N  =  J  d  r  (E  +  B  )/(4h  >>  1,  where  2  y  c  k^  is  the  laser  frequency,  and  h 

is  Planck's  constant. 

This  paper  presents  an  analysis  of  the  transition  from  the  incoherent 

2  Q 

radiation  to  the  coherent  radiation  in  an  FEE  oscillator.  The  model  of  the  FEE 
oscillator  is  described  in  Sec.  II.  The  equations  governing  the  complex  amplitude 
of  the  radiation  in  terms  of  the  particle  trajectories  are  derived  in  Sec.  Ill; 
and  the  equations  governing  the  particle  trajectories  in  terms  of  the  radiation 
field  are  derived  in  Sec.  IV.  The  results  of  Secs.  Ill  and  IV  are  combined  in 
Sec.  V  to  obtain  the  self-contained  radiation  dynamics  equations.  The  equation 
describing  the  dynamics  of  the  radiation  energy  rate  matrix  is  derived  in  Sec. 

VI.  The  solution  of  the  energy  rate  equation  is  obtained  in  Sec.  VII.  The  three- 
dimensional  effects  of  the  spontaneous  radiation  are  incorporated  into  the  one- 
dimensional  model  through  a  filling  factor  in  Sec.  VIII.  The  analysis  of  the  FEE 
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oscillator  start-up  process  is  now  completed-  We  examine  a  limiting  case,  in  Sec. 
IX,  where  the  electron  pulse  length  is  long.  In  the  final  section,  Sec.  X,  we 
compare  our  numerical  results  with  Stanford's  KKL  oscillator  data.  Finally  in  this 
section  a  number  of  possible  methods  are  suggested  to  shorten  the  FKL  oscillator's 
start-up  time. 
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1 1 .  KEL  Oscillator  Start  Up  Model 


The  schematic  representation  of  the  FEL  oscillator  model  used  in  our  analysis 
is  shown  in  Fig.  1.  The  resonator  defined  by  plane  reflectors  at  z  =  0  and  z  =  L 
contains  the  wiggler  magnetic  field  located  between  z  =  LQ  and  z  =  LQ  +  Lw-  The 
total  resonator  losses  are  modeled  heur istically  by  a  Q  factor.  The  highly 
relativistic  pulsed  electron  beam  enters  the  resonator  from  the  left  at  z  =  0  with 
axial  velocity  v^  e^.  Within  the  wiggler  field  the  axial  pulse  velocity  is 
reduced  slightly  to  v  e^.  The  electron  beam  pulses  are  spatially  periodic  with 
period  L^.  Though  is  arbitrary  in  the  analysis,  it  is  clear  that  for  proper 
matching  between  the  beam  and  radiation  pulses  that,  should  be  approximately  an 
integer  times  2vQL/c.  The  radiation  pulse  in  the  wiggler  field,  when  overlapping 
with  the  beam  pulse,  can  travel  at  a  velocity  slightly  less  than  c,  and  the  effect 
is  called  laser  lethargy.  It,  therefore,  becomes  necessary  to  slighty  mistune 
(shorten)  the  resonator  length  to  optimize  the  interaction.  This  effect  is  fully 
taken  into  account  and  is  discussed  in  detail  later.  The  axial  profile  of  the 
electron  beam  pulses  are  left  arbitrary  but  have  a  characteristic  length 
2^  <<  L^.  The  entering  electron  beam  is  monoenerget ic  with  no  spread  in  either 
the  longitudinal  or  transverse  velocities.  The  radiation  pulse  is  assumed  to 
undergo  little  change  in  phase  and  amplitude  during  a  single  pass  through  the 
resonator,  i.e.,  low  gain  operating  regime.  The  wiggler  parameters  are  taken  to 
be  fixed  and  space  charge  effects  neglected.  Finally  the  analysis  is  performed  in 
the  small  signal  regime,  i.e.,  to  first  order  in  the  radiation  field. 
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III.  Reduced  Wave  Equation 


We  will  represent  the  radiation  field  within  the  resonator  by  a  superposition 
of  spatial  modes,  which  are  such  that  the  tangential  electric  field  vanishes  on 
the  mirrors.  The  vector  potential  of  the  radiation  field  is  written  as 

CO  i(0  t  A 

A  (z,t)  =  V  a  (t)  sin  (k  z)  e  "  e  +  c.c.  (1) 

~R  •  n  n  x 

n-  l 

where  k  =  (o  /c  =  im/L,  a„(t)  is  the  Fourier  coefficient  of  the  nth  mode,  and  c.c. 
n  n  n 

denotes  the  complex  conjugate.  The  vector  potential  of  the  linearly  polarized 
wiggler  field  is  non-zero  only  in  the  interval  Lq  <  z  <  +  Lw  and  is  taken  to  be 


A  (z)  =  A  cos  (k  z)  e 

~w  W  W  X 


where  k  =  2 rr /l  ,  i  Is  the  wiggler  wavelength  and  |A  I  »  Jar  )  -  The  one- 
w  w  w  x 

dimensional  wave  equation  for  AR,  including  a  phenomenological  loss  term,  is 


^  (Z’° 


(z,t) 


(3) 


where  the  current  density  J  will  eventually  be  taken  to  be  linear  in 

Ar,  v  =  ui^/Q,  uij  is  the  characteristic  laser  frequency,  and  Q  is  the  quality 

factor  associated  with  the  resonator.  In  the  FEL  the  characteristic  laser 

2—1/2 

frequency  is  =  (1  +  PQZ)  Yoz  voz  kw>  where  ?oz=  voz/c  and  Yoz=  (1  -  l\,2) 

In  (3)  the  Q  is  defined  in  the  usual  way  such  that  in  the  absence  of  a  driving 

2 

current  the  electromagnetic  stored  energy  (proportional  to  la^COT)  decays  like 
exp  (-u),  t/Q).  Note  that  in  (3)  it  is  assumed  that  all  the  significantly  excited 
longitudinal  modes  have  the  same  Q. 

The  actual  discrete  beam  density  is 
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n(z,t)  =  -^  \  5  [z  -  z  (z^,  O]  (4) 

and  o,  is  the  cross  sectional  area  of  the  electron  beam  and  z  (z  ,  t) 
b  o  j 

represents  the  axial  orbit  of  the  jth  electron.  At  t=0  the  initial  axial  position 

of  the  ith  electron  is  z„,,  i.e.,  z  (z  t=0)  =  z  .. 

°J  oj’  oj 

The  fluid  like  beam  density  can  be  defined  as 

nQ(z,t)  =  <  n(z , t )  >  (5) 

where  the  bracket  <  >  denotes  the  ensemble  average  of  the  enclosed  quantity. 

The  ensemble  average  in  (5)  is  over  uncorrelated  charged  sheets  (electrons). 

Using  (4)  we  note  that  the  ensemble  average  of  the  density  n(z,t)  is 

<  ~T  I  S(z  "  z  (z  ,.0  >  =  /  dz  n  (z  ,0)  S(z  -  z(z  ,t))  -■*  n  (z,t)  (6) 

Oi  >  i  '■'J  uuo  u  o 

b  j-l 

where  n^z^.O)  Is  the  initial  spatial  density  distribution  of  particles.  The 
fluctuating  part  of  the  density  is  given  by  n(z,t)  -  <n(z,t)>.  The  effective  non¬ 
linear  driving  current  density  is  given  by 

J  (z,t)  =  J^Cz.t)  +  Jinc(z,t)  (7a) 

where  is  the  cohgrent  current  driving  the  stimulated  radiation  (gain) 

and  J.  is  the  incoherent  contribution  due  to  the  discrete  nature  of  the 
~inc 

electrons  and  is  responsible  for  the  spontaneous  radiation  (shot  noise).  The 
coherent  and  incoherent  current  densities  are  respectively  given  by 

J  =  -  lelv  F  <n( z , t ) >  (7b) 

~c  ~w  c 
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and 


J.  =  -  I e I v  F  [ n( z , t )  -  <n(z,t)>] 
'-inc  'W  Inc  1 


(7c) 


where 


v  =  c  B  =  lei  A  ( z  )  /  y  m  e 
— w  -mw  ~w  o  o 


v  cos  (k  z)  e 

W  W  X 


(8) 


Is  the  wiggle  velocity  defined  over  the  region  L  <  z  <  L  +  L  , 

o--ow 

2  2-1/2 

v  =  I  e  |  A  /(y  m  c)  and  y  =  (1  -  v  ~/c")  •  The  usual  filling  factor  associated 

w  w  'o  o  o  o  b 

with  the  coherent  radiation  is  Fc  =  o^/o  »  where  o  is  the  transverse  area  of  tlu 

resonator  radiation  mode.  The  filling  factor  associated  with  the  incoherent 

radiation  is  written  as  F .  =  / f  /  f .  The  term  f  is  a  loss  factor  due  to  the 

inc  m  m 

a 

finite  size  of  the  mirror  at  z=L  and  is  given  by  f  =  [  2y  r  /(I  +  y  B  )L1“ 

m  'o  m  'o  w 

where  r  is  the  mirror  radius.  In  obtaining  f  we  have  taken  the  incoherent 
m  m 

radiation  divergence  angle  to  be  =  ( 1 / yq  +  Bw)«  The  origin  of  the  second  term  in 

the  expression  for  F|nc  arises  from  the  one  dimensional  statistics  performed  on 

the  uncorrelated  particles.  In  Sec.  VIII  we  show  that  this  term  is  given  by 
6  2  -1-2 

7.  =  Ou(y  /y  )  /  ( X.  Y  )~  wtiere  X.  =  l  ( 1  +  6  )  y  is  the  characteristic  laser 

b  o  oz  l.oz  L  w  oz  oz 

wavelength. 

It  should  be  emphasized  that  in  our  model  the  electrons  are  actually 

represented  by  sheets  of  charge.  The  surface  charge  of  each  sheet  is  -|el/o,  and 

b 

the  sheets  (electrons)  are  taken  to  be  uncorrelated. 

To  obtain  an  equation  for  an(t),  the  Fourier  coefficients  of  the  radiation 

field,  we  first  substitute  (1)  together  with  (7)  and  (8)  into  the  wave  equation 

(3).  Taking  an(t)  to  he  a  slowly  varying  function  of  time, 

t.e.,  la  /a  I  <<  ,jj  ,  there  results  on  neglecting  small  terms 
n  n  n 


00  a)  ioj  t 

-2i  Y  — —  fa  (t)  +  ^  a  (t )  1  sink  z  e  n  +c.c. 

2  L  n  2  n  1  n 

n=L  c 

F  oo 

=  4-n  |  e  I  6  cop  k  z  f  ■  Y  6(z-z(z  ,,t))  -  (F.  -  F  )  n  (z,t)] 

w  w  k  a,  .s  or  inc  c  o  ‘ 

b  j-1 


0(z  -  L  )  0  (L  +  L 
o  o  w 


z) 


(9) 


where  g  =  v  /c  is  the  normalized  wiggle  velocity,  0(x)  is  the  usual  Heaviside 
w  w 

unit  step  function,  and  the  dot  denotes  a  time  derivative.  By  multiplying  both 
sides  of  (9)  by  sin  kffl  z,  integrating  over  z  from  0  to  L,  and  keeping  the 
appropriate  resonant  terms  we  obtain 


irlel  v  L  i(k  +  k  )z  -iw  t 

♦..-V.,,  w  ,  ,  n  w  n 

3n(t)  =  —  an(°  +  Tl -  /  dz  6 

0 


F. 

r  inc 


j  =  l 


l  5(2  "  z^zoj,,:))  "  (Finc*  Fc>  no(z*°) 


9(z-L  )  0(L  +  L  -  z ) . 
o  o  w 


(10) 


To  evaluate  a^Ct),  knowledge  of  the  axial  orbit,  i.e., 


z(z 


is  required. 
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IV.  Particle  Dynamics 


9 


z(  z 


Z  .  +  V  t 

O  J  o 


(1A) 


where  t  <  (L  -  z  ,)/v  .  Within  the  wiggler  the  trajectory  of  the  jth  electron  is 
o  oj  o 


zUo-.t)  =  z(o)(zoj,t)  +  5?(zoj,t) 


(15) 


where  z^°^(z  ,,t)  =  v  z  ,/v  +  (l  -  v  /v  )  L  +  v  t  is  the  unperturbed  orbit 

oj  oz  oj  o  oz  o  o  oz 

and  5 z  is  the  displacement  due  to  the  ponderonot i ve  force.  Equations  (12)  and 

(15)  are  valid  for  times  such  that  the  particle  is  in  the  wiggler,  i.e., 

(L  -  z  ,)/v  <  t  <  (I.  -  z  ,)/v  +  L  /v  .  Substituting  (15)  into  (12)  and 

o  o j  o  -  -  o  o  j  o  w  oz 

linearizing  we  find  that  the  longitudinal  displacement  of  the  jth  electron 
sat  isf ies 


'z  (z  .,t) 
OJ 


I  e  I  B  It 

w  w 

2v  m 
o  o 


■>  r 


(t)  e 


-i  (k  +  k  ) (z 
n  w 


'  .+  L  -  L" ) 

0  J  o  o 


-ipnt 

e 


+  c .c  .  !0 


z  I.  I.  -  z  .+  L' 

_°i _ 2_  +  t)  o  f~ _ ^ _ « 


t) 


(16) 


where  ft  =  v  /c.  u  =v  (k  +  k  )  -  u  =  v  k  -  c  k  (1  -  p  )is  the  frequency 
oz  oz  n  oz  n  w  n  oz  w  n  oz 

mismatch,  z'  .  =  v  z  ,/v  ,  L'  =  v  h  /v  and  L'  -  v  L  /v  .  We  now  invoke  the 
oj  ozojo  o  ozoo  w  o  w  oz 

low  gain  assumption  by  taking  the  coefficients  nn(t)  to  be  constant  during  the 

tine  the  jth  electron  is  within  the  wiggler  region.  Integrating  (16)  twice,  using 

the  low  gain  assumption,  and  taking  the  initial  conditions  such  that  the  relative 

displacement  and  relative  displacement  velocity  are  zero  at  the  entrance  to  the 

wiggler,  i.e.,  5z  =  Bz  =  0  at  t  =  (L  -  z  ,)/v  ,  we  find  that 
"  o  oj  o 


^(Zo.,t)  =  f 


lei R  k  ®  ,  -i (k  +  k  )(z' .+  L  -  L') 

w  w  ,  i  ,  ,  -2  m  w  oi  o  o 

a -  1  N  e 

'Y  m  .mm 

o  o  m=  1 


10 


)  -l) 


1  +  c.c.  (17) 


-iut 

m 


+  (iP  (t 


(Lo" 


*0^ 


-ip  (L  -  z  . ) / v 
mo  oj  o 


where  expression  (17)  is  valid  for  times  such  that  (L  -  z  ,)/v  <  t 

o  cj  o  - 

<  (Lq  +  -  z^j  ]/vq  anc^  is  zero  prior  to  this  time  Interval.  Expression  (10) 

together  with  (17)  describes  the  linear,  low  gain,  longitudinal  dynamics  of  the 
jth  particle  within  the  wiggler  field. 
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V .  Radiation  Dynamics 

We  now  return  to  the  evolution  of  the  radiation  field.  Substituting  (15) 
together  with  (17)  into  (10),  introducing  coefficients  ^(t)  =  k^a^t),  and 
expanding  the  delta  functions,  the  expression  for  the  time  rate  of  change  of  the 
Fourier  coefficients  is  given  by 

6  (t)  =  -  %  b  (t)  +  S  (t)  +  R  (t)  (18) 

a  2  n  n  n 


where 


S  (t) 
n 


it  I  e  I  v 
_ w 

L  °b 


F 


inc 


L  +  L 

r  O  W 


L 

o 


dz 


i(k  +  k  )z  -  iu)  t 
n  w  n 

e 


(  I 

J-l 


6(z-z^°^ (ZQ j , t) )  -  obnQ(z,t)), 


(19a) 


Rn(t)  = 


it  I  e  I  v  F 
w  c 

L 


L  +  L 
o  w 


L 

o 


dz 


l  (k  +  k  )z  -  iu)  t 
n  w  n 

e 


n 

o 


(z,  t  ) 


(19b) 


and  6z(zoj,t)  is  given  by  (17).  On  the  right  hand  side  of  (18),  the  first  term 
represents  the  resonator  loses,  S  (t)  represents  the  spontaneous  or  incoherent 
radiation  term,  and  the  stimulated  or  coherent  radiation  is  represented  by 
Rn(t).  Substituting  the  linearized,  low  gain,  longitudinal  orbit  of  the  jth 
particle  within  the  wiggler  field  given  by  (17)  into  (6)  the  stimulated  terra  in 
(18)  can  be  expressed  as 


Rn(t) 


00 


V 

t. 

m=l 


G 

nm 


(t)  b  (t) 
m 


(20) 
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where 


G 

ntn 


(t) 


l  n  le  1 2v  2k  F  i  ( u  “  u  )t  i(k  -k  )(L  -  L  ) 

w  c  n  m  nmoo-z 

—z - ; -  e  e  u 

2y  ra  cl  m 

'o  o 


L  +  L'  -v  t 


°  /  W  °  dz  n  (z  ,0) 
J  o  o  o 


L (k  -  k  )v  z  / v 
n  m  oz  o  o 


L  -v  t 
o  o 


i1  +  'l^(t 


lvJm^ 

(V  Zo)/vo)  “  11  e 


(L  - 
o 


z 

o 


)/v  ] 

0  \ 


(21) 


where  L  =  L  v  /v  L  •  The  time  rate  of  change  of  the  Fourier  amplitude  given 

w  w  o  oz  w 

in  (18)  can  therefore  be  put  into  the  form 


b  (t)  =  S  (t)  + 
n  n 


00 


l 

m=l 


(G  (t)  -  %  S  ) 
'  nm  2  nm; 


bm(c) 

m 


(22) 


where  S  (t)  is  the  spontaneous  radiation  source  term,  G  (t)  bm(t)  represents  the 
n  nm  m 

dielectric  response  or  gain,  6  is  the  Kronecker  delta,  and  (v/2)  6  b  (t)  is 

nm  nm  m 

the  loss  term  due  to  the  finite  Q  of  the  resonator.  The  matrix  G  defined  by  the 

elements  G  (t)  will  be  referred  to  as  the  gain  matrix, 
nm 
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VI .  Derivation  of  Energy  Rate  Equation 

The  total  ensemble  average  electromagnetic  energy  within  the  resonator  is 


W 

em 


(t) 


/  d3r  <  K2  +  B 2 > / 8 ti 
vol 


oRl 

zr 


V  <  h  b  *  > 

n  n 


(23a) 


and  the  electromagnetic  power  flowing  axially  within  the  resonator  is 


P 

'em 


(z,t) 


c 


4  TT 


/  dA  <  E  x  B  > 
area 


co. 


1  6n 


00 

\ 

n  ,m=l 


<  b  b  *  > 
n  m 


-i(k  -  k  )  (z  -  ct) 
r  n  m 


i(k  -  k  )  (z  +  ct) 
n  m 
e 


e  +  c  .  c  . 

z 


(23b) 


where  the  bracket  <  >  denotes  the  ensemble  average  over  nncorrelated  sheets 

(electrons)  of  the  enclosed  quantity.  From  (23a)  and  (23b)  it  is  clear  that  the 
quantity  of  real  interest  is  the  energy  density  matrix  £  defined  by  the  elements 


e  (t)  =  <bn(t)  b 
nm  n  m 


(t)>. 


(24) 


In  terms  of  the  energy  density  matrix  in  (24),  the  total  electromagnetic  energy, 
Wem,  and  the  electromagnet ic  power,  P^  (z,t),  are 


W  (t) 
em 


°rl 


V 

n=l 


(t) 


(2Sa) 


and 


P 

~em 


(z,t) 


1 6n 


oo 

\ 


n,m=i 


e  (t) 
nm 


-  k  )  (  z  -  ct) 
m 


i(k  -  k  ) (z  +  ct ) 
n  m 
e 


e  +  c  .  c  . 
z 


(25b) 
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We  now  derive  the  rate  equation  for  the  energy  density  matrix.  Writing  (22)  In 
vector  notation  yields 


i  (t)  =  S  (t)  +  [C  (t)  -  j  I]  b(t ) 


where  I  is  the  unit  matrix.  Solving  (26),  with  initial  condition 
we  obtain 


b  (t)  =  /  X  (t)  X  l(t')  S(t')  dt 

0 


where  X  (t)  is  defined  by  the  equation 


*  (t)  =  [G(t )  -  |  I  1  X  (t) 


with  initial  conditions  X(0)  =  l.  The  energy  density  matrix  is 

£  (t)  =  <  b(t)  b"(t)  > 

-l  3  2  2, 

where  Trace  (g)  =  (oRh)  J  d  r  <  E  +  B  >/2  and  the  superscript 
Hermltlan  conjugate.  Using  (27)  together  with  (29)  we  find  that 
the  rate  equation 

g  (t)  =  fg  (t)  -  v  I]  g  (t)  +  Y  (t)  +  11. C. 

where  )’  (t)  =  /  <  S  (t)  if  (t')  >  (X  (t)  X_l(t'))H  dt' 

»  0 

and  H.C.  denotes  the  Hermltlan  conjugate  of  the  preceding  terms. 


(26) 

b(0)  =  0, 

(27) 


(28) 


(29) 

H  denotes  the 
■  (t)  satisfies 

(30) 


It  can  be  shown 


15 


that  the  ensemble  average  of  S,,(t)  S  (t)  can  he  expressed  as 

n  m 


it  le  I 


V  S  (t)  S 
n  m 


(t')  >  = 


2  2 

F 7  ill  t  -tg  t' 
v  Inc  n 
-  e  e 


i(k  -  k  )(L  -  L') 
n  m  o  o 


L  +  L'  -v  t 
o  w  o 

f  dz  n  ( z  , 0 )  e 
o  o  o 

I,  -v  t. 
o  o 


i(k-k)v  z/v 
n  m  oz  o  o 


Il(z  -  L  +  v  t')  0(L  +■  1.  '-z  -  V  t') 

o  o  o  o  u  o  o 


(31) 


By  noting  the  limits  of  integration  as  well  as  the  arguments  of  the  Heaviside 
functions  in  (31),  it  is  clear  that  the  ensemble  average  <  S^(t)  Sm(t)  >  is  non¬ 
zero  only  for  t  -  t "  <  L  / v  ,  i.e.,  when  there  is  an  electron  sheet  in  the 

w  oz 

wiggler.  It  can  he  shown  that  in  this  interval  X  (t)  X  1  (t')  =  I  if  tie  gain  per 
pass  is  somewhat  less  than  unity.  Hence  the  source  term  of  the  energy  rate  Hq. 
(30)  is  simplified  to 

t  ~ 

\  (t )  »  /  <  S(t)  S  (t')  >  dt".  (32) 

=  0 


The  elements  of  (32)  take  the  form 


it  lei  v  2  K  i  ( n  —  u  )  t  i(k  -  k  )(L  -  L') 

V  ( t )  =  if  _ -  )  liu.  e  n  ,n  e  11  m  °  ° 


(33) 


L  +  l'  -v  t 
o  w  o 


i  ( k  -  k  )  v  z/v  ,  i n  ( t  -  ( L  -  z  )/v) 

/  .  ,  n  m  oz  o  o  -1  ,  in  o  o  o 

J  dz  n  (z  ,0)  e  1 1  1  -  e 

'  o  o  o  m 


L  -v  t 
o  o 


This  completes  our  formal  derivation  of  the  energy  rate  equation  given  by  (30). 
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Due  to  the  complicated  structure  of  both  the  sain  matrix  (21)  as  well  as  the 


spontaneous  source  matrix  (33),  it  is  convenient  to  further  reduce  these  terms  to 
a  more  manageable  form.  To  this  end  we  define  a  time  variable  t^,  such  that  tN  is 
the  time  that  the  center  of  the  Nth  electron  pulse  enters  the  wiggler  field,  that 
is 


tN  =  ((N  -  !)  hb  +  Lo)/vo  (34) 

where  N  is  a  non  zero  positive  integer.  During  the  electron  pulse  propagation 
through  the  wiggler  field,  the  independent  time  variable  is 


t  =  tN  +  T 


where  0  <  r  <  L  /v  .  To  simplify  the  gain  matrix  and  spontaneous  source  matrix 
-  w  oz 

in  (21)  and  (33)  we  note  that  these  matrices  involve  integrals  of  the  general  form 


l  (t)  = 
nm 


L  +  l."-v  t  .  ,  , 

o  wo  t (k  -  k  ) v  z  / v 

|  dz  n  (z  )  e  1  (t,z  ).  (35) 

,  1  o  o  o  nn  o 

t.  -v  t 
o  o 


The  generic  integral  in  (35)  can  be  evaluated  for  two  representative  electron 
pulse  shapes  of  characteristic  width  4^  given 


I 

nm 


<v 


t) 


(N-l)LhJno4b 


-i(kn- 


k  )  (N 
m 


l)v  L  / v 
oz  b  o 


Pn_  (3b) 
nm 


where 


’’nm 


/  Tt 
2 


-((k  -  k  )  9/U)2 
v  n  m  b  '  ^  r , 

e  ,  Gaussian  profile 


sln( (k  -  k  )  Z./2) 
n  m  d 


(kn"  km) 

n  m 


'Vr 


(37a, k) 


,  square  profile. 


The  expression  in  (37a)  is  for  a  Gaussian  electron  beam  pulse  shape,  i.e.. 


n  (z  )  =  n  e 
o  o  o 


■(2  z  /O' 

O  D 


(38a) 


while  the  expresssion  in  (37b)  is  for  a  square  pulse  shape,  i.e., 


n  (z  )  = 
o  o 


n  ,  /2  <  z  <  i  /2 

o  b  -  o  -  b  , 

0,  otherwise. 


(38b) 


Using  the  result  contained  in  (36),  together  witli  (34),  both  the  gain  matrix  and 
the  spontaneous  source  matrix  can  be  reduced  to 


.  1.  u>r  ,  -i(k  -k  ) (N— 1 )  L. 

»  *  1  b  „  2  f  n  m  b 

C  (tM+  t)  =  q-  t—  — —  0  k  c  F  c 

nmN  8L  vwwc 

o 


u  p  i(u  “  u  )  r 
nm  nm  n  m 


iu  f 

(1  +  (in  x  -  1)  e  1 
v  m 


(39a) 


and 


\  ft  +  i  =  — 

■  nm  \  T  4  L,  La,  "w  ‘  inc 


uv"  m  -  ,  -i(k  -  k  )  (M-l)  L. 

o  o  2  „2  vn  m  '  b 

v  F 


«  P 
nm  nm 


i(u  ”  u  )  t  iu  t 

n  m  , ,  m 

e  ( 1  -e 


(39b) 


where 


a  =  exp  (~i(k  -k  )(l-p  )6  *L  ) 
nm  n  m  o  o  o 


2  2 

and  =  4u|e|  n^/m^  is  the  peak  beam  plasma  frequency.  In  obtaining  (39a)  and 

(39b)  we  replaced  v  l./v  by  t,. 

oz  b  o  1  h 
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In  the  absense  of  "laser  lethargy"  exact  resonance  between  the  electron  beam 
pulses  and  the  radiation  pulses  occur  when  the  mirror  separation  is  equal  to 
Lb/(23Q)  where  8q  is  the  normalized  axial  pulse  velocity  outside  the  wiggler 
field.  This  condition  implies  that  the  round  trip  of  the  radiation  pulse,  if  it 
were  traveling  at  c,  equals  the  beam  pulse  period.  However,  since  the  radiation 
pulse  velocity,  in  the  wiggler  region  when  overlapping  with  the  electron  pulse,  is 
slightly  less  than  c,  it  is  necessary  to  have  the  mirror  separation  slightly  less 
than  (L^/2p^)  for  optimum  overlap  of  the  beam  and  radiation  pulses.'’  ^  With  this 
in  mind  we  define  the  mirror  separation  to  be 


L  =  L  +  &L 
m 


(40) 


where  L  -  L  /(2f?  )  >>  I SL I .  In  (38)  and  (39)  the  only  term  sensitive  to  slight 
mho 

variations  in  the  mirror  separation  is  the  common  leading  term  exp  ( — i ( k  -  k  ) 

n  m 

(N-l)Lb).  Substituting  (40)  into  (39a)  and  (39b)  and  assuming  SL  small  we  find 
that  the  gain  and  source  matrix  elements  become 


S  (t,  ^ 

nra  N  8  L  y 


2 in  (n-m )  ( N  —  1 )  .S L/L 


k  c  F  a  p  e 
w  c  ran  nm 


i(u  -  m  ) t 
e  n  m 


1  M  t 


1  f  (i|i  i  '  1)  r 
‘  m 


(41) 


and 


9  .  ,9 


V  7T  g  n  v 

\  _  'b  o  w  2  2i  tt  (n-m)  (N-l )  (SL/L 

’  (t.T+  t)  *  —  - : -  a  o  F  e  m 

•nm  N  L  I,  a.  nm  nm  i  no 


T  G 


i(u  “  li  )t  t 

n  m  Ip  t/Z 


sin  \i  t  /  2 
n 

”ii  t/2 
n 


f- 


-1. 


(42) 
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The  rate  of  change  of  the  field  energy  density  matrix  given  in  (30),  together 
with  the  expressions  for  Gnm  and  ^  in  (41)  and  (42),  can  be  still  further 
reduced  by  invoking  the  low  gain  per  pass  approximation.  The  low  gain  per  pass 
assumption  implies  that  §  changes  slightly  during  a  single  pass  of  the  radiation 
pulse.  Hence,  by  taking  fi(tN+  t),  where  r  <  VVoZ’  to  be  nearly  equal  to  fi  <£„> 
on  the  right  hand  side  of  (30),  we  can  integrate  (30)  together  with  (41)  and 
(42).  Doing  this  we  find  that  the  elements  of  g  at  time  tN+  t  are  given 
approximately  by 


I)  -  <i  -  »  <>  e„„«s)  +  Vr'1 


♦  ’  lC„!  <tS’,)  W  +  ent(V  C 


,s3) 


il=l 


where 


li  —  —  B  2  k  c  F  e 

32  L  yo  w  w  c 


2in  (n-m)  (N-l  )  6L/I. 


m  a  p  g  (t)  (44a) 
nm  nm  nm 


and 


tt2  I e  [  2n  v  2  ,  2in(n-m)(N-l)«L/L 

0  W  F  ^  e  no 

2L  - WT  inc  nm  nm  nm 


*Tj 


h  (t)  (44b) 


where 
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g  (  r )  = 
nm 


3  ix  x 

r  (  n  f .  .  m. 

- 7T  c  ( l  +  — )  sin  x 

2  '  x  n 

xx  n 

n  m 


t(x-x)  s in( x  -  x  ) 

n  m  ,  n  m  i 

.  x  - -  - 

n  x  -  x 
n  m 


x  e 
m 


and 


(43a) 


T*  2  2  2 

h  ( t )  —  -  i  sin  x  +  s i n  x  -sin  (x-x) 

nm  xx'  n  m  n  m 

n  m 


(43b) 


-i  sin  x  cos  x  -  sin  x  cos  x  -  sin  (x  -  x  )  oos(x  -x  ) 
n  n  n  m  n  m  n  n 


and  x  =  u  t/2  =  Iv  k  -  ck  (1  -  v  /c)lr/2.  Note  tliat  since  is  ilermilian,  so 

n  n  oz  w  n  oz  1  nm 

is  the  spontaneous  source  matrix  Snn  in  (44b).  Tlie  fact  that  Snn  is  Hermitian  is 

simply  a  consequence  of  ttie  fact  that  ■  by  definition  is  Hermitian,  (see 

(24)).  Setting  r  =  l.  / v  in  (43)  gives  t lie  energy  density  matrix  after  the  Nth 

w  oz 

beam  pulse  has  tranversed  the  wiggler.  The  results  obtained  by  numerically 
solving  (43)  for  various  experimental  parameters  will  be  presented  later. 
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VIII.  Spontaneous  Radiation  Source  Term 

The  spontaneous  radiation  source  term  in  (44b)  has  been  obtained  from  a  one 
dimensional  analysis  of  the  wave  equation.  Because  of  the  one-dimensional 
character  of  the  analysis  the  spontaneous  source  term  does  not  properly  represent 
the  incoherent  radiation  source.  A  proper  three-dimensional  treatment  of  the 
spontaneous  radiation  is  necessary  to  properly  consider  the  statistics  of  discrete 

uncorrelated  particles  as  well  as  to  separate  the  "velocity"  and  "acceleration" 

30 

(radiation)  electromagnetic  fields.  The  present  one  dimensional  treatment 

represents  the  electrons  as  uncorrelated  charged  sheets  and  not  as  point 

particles.  To  correct  for  the  one-dimensional  limitations  of  our  analysis  of  the 

spontaneous  source  term  we  have  included  in  the  incoherent  current  density  (7c),  a 

filling  factor  which  contains  the  term  /'.  .  This  term  is  included  so  that  the 

total  emitted  spontaneous  radiation  agrees  with  the  well  known  value  obtained  from 

Lartnor's  formula.  We  have  justified  this  procedure  by  performing  a  proper  three- 

dimensional  treatment  of  the  spontaneous  source  term;  this  three-dimensional 

analysis  will  be  published  elsewhere.  To  obtain  the  factor  £  in  the  spontaneous 

source  term,  i.e.,  in  the  filling  factor  F.  ,  we  compare  the  total  emitted 

l  nc 

radiation  energy  from  (43)  with  that  obtained  from  Larmor's  formula  with  the  loss 
terms  f  set  equal  to  unity  and  v  =  0.  From  (43),  the  diagonal  elements 
of  c  satisfy 


enn(tN  +  0  =  2  ^ 


2,  ,2  2 
n  lei  n  v 
o  w 


s  <  n  u  T 
n 


b  o. 


(46) 


where  we  have  used  the  expression  for  '  in  (44b)  and  are  considering  a  square 

■  nm 

shaped  electron  beam  pulse,  i.e.,  p  =  1.  We  now  want  to  compare  (46)  with 
I, armor's  radiative  formula.  The  total  instantaneous  power  radiated^  from  a 
single  particle  is 
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pL  2  lei “  f,-2  .  ... 

P  "a  — T~  ~  C  v  x  v)). 

em  i  i  o  ~  ~  ~ 

c 


The  velocity  of  a  single  particle  in  the  uiggler  is  v  *  ^  oz  e^  +  v^  cos 

(k^v  t)  e^.  Using  (47)  we  find  that  the  total  energy  radiated  during  a  time 

t  <  L  /v  ,  by  a  beam  pulse  consisting  of  £,  n  o,  particles,  is 
-  w  oz  bob 


2 

,,L  l  lei  Yo  .  .  .2 

em  =  3  - - —  Vo°b  (3oz  Vw)  T> 

^oz 


The  total  spontaneous  electromagnetic  energy  within  the  resonator  is  given  by 
(25a)  with  oR  replaced  by  o^.  Substituting  (46)  into  (25a)  gives 


a,  L  » 


<  L  /v 


W  (t,,  +  T)  = 
etn  N  4tj 


f  dn  f  dT  W  °Z 


2,2  2 
2 1,  it  I  e  I  n  v 


sin  ii  t 


V  ^  A.  II  U  l 

o  w  /l  x  n  ,  ,,  ,  . 

-  F,  - —  +  W  (t.,) 

inc  u  em  N 

n 


where  we  have  approximated  the  sum  by  an  integral.  Integrating  (49)  over  t  and  n 
gives 


.  ,2  ?  n  v 

em  N  r  em  2c  (l  -  g  )  inc  T’ 

07. 


Comparing  (48)  and  (50)  we  find  that,  for  f  =  1, 

m 


F2  =  r  ,  f  ±  l6  _!> _ 

inc  1  y  ,2 

02  (VozJ 


where  X,  is  the  laser  wavelength,  X  =  i  /  ( 1  *  4  )  y 

L  B  L  w  nz  oz 
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IX.  Long  Beam  Pulse  Limit 


A  limiting  case  which  can  be  fully  evaluated  analytically  Is  that  of  a  long 

pulse  beam,  i.e.,  8.^  <  L.  Though  this  limit  is  not  necessarily  directly 

applicable  to  either  planned  or  completed  pulsed  beam  FKL  oscillator  experiments 

it  does  represent  an  Interesting  limit  of  the  more  realistic  configurations.  If 

the  electron  pulse  widths  are  comparable  but  somewhat  less  than  the  mirror 

separation,  L,  the  gain  matrix  as  well  as  the  spontaneous  source  matrix  in  (44a) 

and  (44b)  approach  a  diagonal  form.  This  can  be  seen  by  noting  that 

for  i.  <  L,  the  matrix  defined  by  p  and  used  in  (44a)  and  (44b) 
b  ~  nm 

approachs  (/n/2)  6  for  a  Gaussian  beam  pulse  and  6  for  a  square  shape  beam 
nm  nm 

pulse  where  6  is  the  Kronecker  delta.  The  diagonal  form  of  (44a)  and  (44b)  is 
nm 

reasonable  in  this  limit,  since  it  is  the  off  diagonal  elements,  in  particular  the 
term  exp  (2ni (n-tn) (N-l ) SL/L^) ,  which  are  responsible  for  the  laser  lethargy  effect 
and  when  the  beam  widtti  is  sufficiently  long  this  effect  is  unimportant.  In  this 
limit  a  single  longitudinal  mode  analysis  would  suffice. 

Therefore,  the  energy  rate  equation  in  (43)  together  with  (44a)  and  (44b), 
for  long  beam  pulses,  takes  the  form 

e  (t  +  x)  »  (1  ~  v  T  +  g  (t))  e  (t  )  +  S  (t)  (52) 

nn  N  nn  nn  N  nn 

where  the  diagonal  gain  and  source  matrix  elements  are  respectively 


g  (t)=G  (tv+  r)  +  G  (t.,+  t) 
nn  nn  N  nn  N 


J.  a  2v  F  3  a  rstn  Xn,2 

16  Ly  ®w  wC  cT  3x  ^  x  ^ 

o  n  n 


Snn(T>  f 

nn  L 


2  2 

iL  it  lei  n  v  ^ 

d  o  w  „2  /  n-v  2 


La, 


sin  x  2 


(53a) 


(53b) 
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and  =  u^t/2.  In  obtaining  (53)  we  have  assumed  a  square  pulse  shape.  Vote 

that  in  g  (t)  and  S  (r),  r  ranges  from  0  to  L  /v  .  Since  c  changes  little 
6nn  nn  *  b  w  oz  nn 

from  pulse  to  pulse  we  may  transform  (52)  into  a  first  order  temporal  differential 

equation.  Since  t.,,,  =  t.,  +  L.  /v  .  (52)  can  be  written  as 
N+L  N  bo 


de  (t) 
nn 


!  g  /At  -  vie  (t)+S  /At 
1  nn  1  nn  nn 


where  At  =  2L,/v  ,  g  and  S  are  to  be  evaluated  at  t  5  1  /v  and  we  have 
b  o  nn  nn  w  oz 

replaced  the  discrete  time  parameter  tjg  with  the  continuous  parameter  t. 
Intregrating  (54)  yields 


c(t) 

nn 


S  (g  -  vAt  )t/At 

nn  ,  nn  , 

-  f  e  -1 


g~  vAt 
nn 


where  enn(t=^)  =  0*  For  times  less  than  a  growth  time,  i.e., 
t  <  At/(gnn~vAt), 


e  (t)  =  S  (t/At  +  (g  -vAt)t2/2At2). 
nn  nn  nn 
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X. 


Numerical  Illustrations,  Experimental  Comparison  and  Discussion 

Our  numerical  Illustrations  are  oirected  towards  a  comparison  of  the  FKL 

oscillator  experimental  results  reported  In  Ref.  4.  In  addition,  we  suggest 

methods,  which  could  substantially  shorten  the  oscillator  start-up  time. 

The  parameters  of  Stanford's  FF.L  oscillator  is  given  in  Table  I.  In  the  FF.L 

oscillator  experiment  a  helical  wiggler  field  was  used.  Since  the  present 

analysis  assumes  a  linear  wiggler  it  becomes  necessary  to  multiply  B  in  Table  I 

w 

by  /2  in  order  to  be  consistent.  The  peak  power  within  the  resonator  as  a 

function  of  the  number  of  beam  pulses  that  have  passed  through  the  resonator  is 

shown  in  Fig.  2  for  six  values  of  the  resonator  mismatch  length  6L  =  L  -L.  / 2 g  . 

m  b  o 

Figure  3  shows  the  asymptotic  gain  as  a  function  of  6L.  The  mirror  mismatch  61.  = 

-3 

-  1.1  x  10  cm  corresponds  to  maximum  gain  but  not  maximum  satruated  power. 

_3 

Maximum  saturated  power  occurs  for  6L  between  0  and  -t.l  x  10  cm.  The  range 

_3 

in  6L  for  nonzero  gain  is  -  3.0  x  10  cm  <  6L  <  0,  in  fair  agreement  with  the 

_3 

experimental  range  of  2.5  x  10  cm.  The  maximum  calculated  multi-mode  (finite 
beam  pulse)  power  gain  is  0.16  whereas  the  single  mode  (continuous  beam)  yields  a 
value  of  0.25.  Finite  beam  pulse  effects  therefore  reduce  the  linear  gain  by 
approximately  60%.  The  maximum  experimental  gain  is  0.10. 

Figure  4  shows  the  spatial  distribution  of  the  electron  pulse  (square)  and 
the  radiation  power  pulses  at  the  entrance  and  exit  of  the  wiggler  for  6L  =  -  1.0 
x  10  cm.  Upon  entering  the  wiggler  the  radiation  pulse  slightly  lags  the  beam 
pulse,  while  exiting  the  wiggler  the  two  are  completely  overlapped.  The 
asymptotic  energy  spectrum  of  the  radiation.  Fig.  5,  is  narrower  and  shifted  with 
respect  to  the  spontaneous  radiation  spectrum. 

F.quation  (43)  suggests  that  one  can  roughly  compute  the  relationship  between 
PN,  the  peak  power  in  the  resonator  after  the  Nth  pulse,  to  P  the  power  emitted 
spontaneously,  by  assuming  a  constant  average  gain  per  pass  g.  An  elementary 
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calculation  yields  when  N  »  1  and  g  <<  1,  Pjj/I’q  **  N  -  1  +  (1  +  g)^  =  N  +  exp 
(gN).  Clearly  when  gN  >>  1  the  result  is  very  sensitive  to  small  changes  in  g  and 
N.  If  one  takes  the  experimental  values  corresponding  to  the  maximum  observed 
final  power  of  PN  =  2.7  x  10^  W  within  the  resonator,  N  =  540  and  the  computed 

-9 

spontaneous  power  of  PQ  =  6.5  x  10  W,  one  finds  that  g  =  0.037.  The 
experimental  value  of  linear  gain  is  0.067.  Tn  view  of  the  sensitivity  to  changes 
in  N  and  g  the  results  are  not  inconsistent.  Moreover  this  effective  value  of  g 
is  smaller  than  the  linear  gain  predicted  by  the  present  model  which  is  reasonable 
since  non-linear  effects  and  initial  beam  thermal  effects  must  lower  the  gain. 
Unfortunately  the  currently  available  data  is  inadequate  to  make  other  detailed 
comparisons  with  this  small-signal  theory. 

Our  analysis  suggests  possible  ways  to  substantially  shorten  the  oscillator 

start-up  time  while  maintaining  high  saturated  power  levels.  The  first  approach 

takes  advantage  of  the  fact  that  the  maximum  linear  gain  and  maximum  saturated 

power  occur  for  different  values  of  6L,  which  we  will  respectively  denote 

by  and  61.,.  By  slightly  increasing  the  frequency  of  the  R.P.  accelerating 

field,  u>  ,  ^ur*n8  t*le  start-up  period,  i.e.,  decreasing  the  beam  pulse 

separation,  the  value  of  6L,  could  be  varied  from  an  initial  value  of  61. ^  to  the 

value  of  61.,,  thus,  decreasing  the  start-up  time  while  maintaining  high  final 

power  levels.  The  required  fractional  increase  in  io  is  1 6L  -  6L_|/I.  «  10 

acc  1  7  b 

for  the  parameters  of  Ref.  3  and  4.  The  same  effect  may  also  be  realized  by 
simply  changing  (increasing)  the  mirror  separation  during  the  start-up  period. 
Another  possible  method  of  decreasing  the  start-up  time  would  be  to  simply 
increase  that  part  of  p£nc  associated  with  mirror  losses,  i.e.,  increase  f^.  This 
could  be  accomplished  by  increasing  the  effective  size  of  the  mirror  located  at  z 
=  L.  The  additional  extension  of  the  mirror  would  necessarily  have  a  different 
curvature.  This  last  approach  should  make  it  possible  to  contain  a  far  larger 
portion  of  the  incoherent  radiation. 
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Table  I:  FEL  Oscillator  Parameters  at  Stanford  University 


Beam  Parmeters 


9 


Beam  Energy,  (y  - 

1 )  ra  c~ 
o 

43  MeV 

Total  Gamma,  y 
'o 

85 

Axial  Gamma,  v 

oz 

69 

Peak  Current,  I 

P 

1.3  A 

Pulse  Width, 

1  mm 

Pulse  Separation, 

25.4  m 

Beam  Radius,  r^ 

0.25  mm 

Wiggler  Parameters 


Wiggler  Wavelength,  9.  3.3  cm 

Wiggler  Amplitude  (helical),  B  =  2n  A  /i  2.3  kG 

WWW 

Wiggler  Length,  Lw  5.3  m 


Resonator  and  Radiation  Parameters 


Resonator  Length,  L 

12.7  m 

Resonator  Losses  (round  trip) 

1.5% 

Radiation  Wavelength,  X^ 

3.3  um 

Spot  Size,  rQ 

0.167  i 

Beam  Filling  Factor,  Fc 

0.01  7 

Incoh.  Rad.  Loss  Factor,  f 

*  m 

0.05 

2 

Rayleigh  Length,  nr^/X^ 

2.71  m 
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INSTANTANEOUS  PEAK  POWER  (WATTS) 


n 


for  Stanford  KF.L  oscillator  experiment  with  various  detuning 
parameters  5L. 


30 


RADIATION  POWER  (RELATIVE  UNITS) 


Fig.  4  Radiation  pulse  power  relative  to  the  spatial  profile  of  the  electron 

pulse  (square)  at  the  entrance  of  the  wiggler  (t  =  t^)  and  exit  of  wiggler 

(t  =  tjg  +  'jw/vzo)>  where  "1  >>  1  denotes  the  electron  pulse  number  for  the 

-3 

Stanford  FKI,  oscillator  experiment  with  <SL  =  -1.0  x  10  cm. 
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Fig.  5 


Asymptotic  energy  spectrum  of  the  radiation  pulse  for  the  Stanford  FEL 

-3 

oscillator  experiment  with  6L  =  -1.0  x  10  cm. 
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